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We have constrained the redshift-distortion parameter (3 = Q. ' 6 /b and the real-space power 
spectrum of the IRAS PSCz survey using a spherical-harmonic redshift-distortion analysis 
combined with a data compression method which is designed to deal with correlated param- 
eters. Our latest result, /3 = 0.4 ± 0.1, strongly rules out /3 = 1. 



1 Introduction 

It has been understood for a long time that redshift surveys are systematically distorted by 
peculiar velocities t3 - so called redshift distortions, which can cause the observed distribution of 
galaxies to become anisotropic. On large scales where linear theory applies, redshift distortions 
are characterised by the parameter ^ f5 = - see Hamilton a for a detailed review. 

In a classic paper, Kaiserlll derived a simple formula for the effect of linear redshift distortions 
for a volume limited survey that subtends a small opening angle at the observer (the "distant 
observer" approximation - all lines of sight are treated as parallel). He showed that power is 
boosted by a factor (1 + /3/i|) 2 , where fi k is the cosine of the angle between wavevector and line 
of sight. Hamilton 00 and Cole, Fisher & Weinberg 1 11 analysed all-sky IRAS surveys using a 
method based on the Kaiser formalism, but were forced to break the survey up into sections 
because of the constraints of the distant observer approximation, losing information about the 
largest (and most reliably linear) scales. 

Rather than fit a square peg into a round hole, Fisher et ala and Heavens & Taylor 
(HT) dropped the plane parallel approximation and used a spherical harmonic decomposition 
to match the spherical nature of the IRAS redshift surveys. HT used spherical Bessel functions 
to decompose the density field radially; using eigenfunctions of the Laplacian retains all the 
advantages of Fourier analysis. HT analysed the IRAS 1.2Jy survey, fitting (3 and the amplitude 
of the power spectrum; Ballinger, Heavens_& Taylor i(BHT) extended this analysis to fit the 
shape of the power spectrum. Tadros et a/ta (T99) applied these techniques to the IRAS PSCz 
survey - section Q includes a review of those results. 

a fl m is the contribution to the density parameter from matter. The bias parameter b depends on galaxy type, 
in this paper it refers to IRAS galaxies. 



2 Spherical Harmonic Formalism 



We will briefly review the formalism for spherical harmonic analysis - see HT and T99 for details. 
The density field of the galaxy distribution p(s) is expanded in terms of spherical harmonics, 
Yt m , and a discrete set of spherical Bessel functions, jj>, 

hmn = Ci n Jd 3 S p(s)w(s)ji (k £n s) Y£ m (0, <j)) , (1) 

where the C£ n are normalization constants and k\ n are discrete wavenumbers. 

These observed coefficients can be related to those of the true underlying density field (6e mn ) 

by: 

hmn = (P0) £mn + E E S nn "W^' + PV$'»') 8 lWn , (2) 

i'm'n' n" 

The transition matrices W, V and S describe the effects of the sky mask, the radial se- 
lection function, the linear redshift space distortion and the small scale distortion correction 
respectively; (po)^ mn is a mean term, non-zero because of partial sky coverage. The transition 
matrices are derived and defined in HT and T99. 

A likelihood approach is used to constrain (3 and the real-space power spectrum: 

- 2 In £[D|/3, P(k)] = ln(det C) + DC^D 7 , (3) 

where C =< DD r > and elements of the data vector D are given by Dg mn = [pe mn — (po)emn]/f> 
where p is the mean number density - Gaussian statistics are assumed. Two different parame- 
terisations of P(k) are used: following HT a fixed shape is assumed and the amplitude is fitted, 
following BHT a stepwise maximum likelihood method is used, allowing the power spectrum to 
assume any shape (within bin discreteness). 

3 Data The IRAS PSCz Survey 

The PSCz surveyEl is a redshift catalogue complete down to 0.6Jy over ~ 83% of the sky with 
a total of ~ 15, 000 redshifts - it is the largest all sky survey in existence. As our method is 
very precise, we minimise systematic errors by using a more conservative flux cut (0.75Jy) and 
sky mask reducing the number of redshifts by a factor of roughly two - see T99. 

4 Previous Results 

Both the fixed amplitude (HT) and stepwise P(k) (BHT) methods were applied to the PSCz 
by T99 for modes with k < 0.13/i Mpc -1 . The first method produced (3 = 0.58 db 0.26 and the 
amplitude of the real space power measured at wavenumber k = O.l/i Mpc -1 of Ao.i = 0.42±0.03 
- see Fig. 2 (dashed contours). Freeing the shape of the power spectrum we find the consistent 
results (3 = 0.47 ± 0.16 (conditional error), and Ao.i = 0.47 ± 0.03 - Fig. 1. T99 also carried 
out extensive tests on simulations, and the methods were found to be reliable. In addition, we 
carried out a suite of tests for systematics effects in the data and found the cut catalogue gave 
consistent results. 

The method was restricted to a small range of wavenumber because the likelihood analysis 
involves the repeated inversion of a large n x n matrix; the time required for this process grows 
as n 3 . More importantly, the matrix rapidly became numerically unstable. The results above do 
not strongly rule out either high (~ 1) or low (~ 0.5) values of /3. It would be nice to overcome 
the matrix problem, extend the k-range and reduce the error bar. 
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Figure 1: The real-space power spectrum of the PSCz redshift survey in dimensionless units (from T99). The 
curve is a CDM model with V — 0.2. The estimated value of j3 is 0.47 ± 0.16 (conditional error). 

5 Data Compression 

It is possible to transform the (length n) data vector D to create a new, smaller dataset, 
while retaining most of the information about the parameters of interested. A new dataset is 
constructed which is a linear combination of the original: 



where D' is a new data vector of length n' < n, with a corresponding transformation for the 
covariance matrix. TTH show that the optimal matrix B which minimises the conditional error 
on a single parameter is made up from eigenvectors of fl 



where the comma refers to a derivative with respect to parameter i. The new covariance matrix 
is smaller and also close to diagonal - hence much more stable. 

TTH suggested extending this to multi-parameter problems by constructing a separate ma- 
trix B for each parameter and then using singular value decomposition to combine the matrices 
efficiently. However, Ballinger§ (see also Taylor et al HH) tested this and found that error el- 
lipses/ellipsoids tended to grow along correlation axes - only conditional errors are constrained. 
Instead it was shown that equation (|B|) could be used to constrain the marginal errors by using 
one or more linear combinations of the original parameters which lie directly along the correlation 
axes. 

Unfortunately solving equation @ involves manipulating uncompressed nxn matrices. This 
usually isn't a timescale problem - it need only be done once - but still suffers from numerical 
difficulties. To avoid this problem, we split the data into several sections and compressed them 
separately, while still retaining full information about correlations between modes in different 
sections. Strictly speaking this is slightly less optimal than compressing the whole dataset in 
one go, but it should make a negligible difference in practice and will not bias the result. 

6 New Results and conclusions 

Fig. 2 shows the old /c max = 0.13/i Mpc -1 result together with the new A; max = 0.2h Mpc -1 
obtained using data compression. The original 4644 modes were compressed down to 2278. The 

This, like the Wiener filter, is one of those marvellous data-handling methods which were invented by as- 
tronomers then sent back in time so that signal processing engineers could use them for the past fifty years. 
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Figure 2: Likelihood contours for the two parameters [3 — Q,q' 6 /b and the amplitude of the power spectrum for the 
PSCz survey. The new data-compression fc max = 0.2h Mpc -1 results (solid contours) are plotted along with the 
original fc max = 0.13ft Mpc" 1 results from T99 (dashed contours). Contours are plotted at intervals of 8\n£. = 0.5. 



new, lower, value of (3 = 0.4 ±0.1 is consistent with the previous result but the error ellipse 
is considerably smaller; the best-fit amplitude of the power spectrum is essentially unchanged. 
(3 = 1, corresponding to f2 m = 1, b = 1, is now strongly disfavoured. More details will be in 
Ballinger et al and Taylor et al (in preparation). 

The low value of (3 is consistent with currently popular cosmological models with a low value 
of ^matter if the IRAS bias parameter is close to unity. The value is somewhat lower than that 
from other analyses of the PSCz survey §@@@, but not inconsistent. It is consistent with the 
recent velocity-velocity comparison results from peculiar velocity catalogues E3, but somewhat 
lower than the corresponding density- density value 0. See Willick El for a discussion. 
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